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Abstract 

We study properties of magnetic field-induced Bose-Einstein condensate of triplons as a function 
of temperature and the field within the Hartree-Fock-Bogoliubov approach including the anomalous 
density. We show that the magnetization is continuous across the transition, in agreement with 
the experiment. In sufficiently strong fields the condensate becomes unstable due to triplon-triplon 
repulsion. As a result, the system is characterized by two critical magnetic fields: one producing 
the condensate and the other destroying it. We show that nonparabolic triplon dispersion arising 
due to the gapped bare spectrum and the crystal structure has a strong infiuence on the phase 
diagram. 

PACS numbers: 75.45. +j, 03.75.Hh 
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Bose-Einstein condensation (BEC), a macroscopic quantum phenomenon, occurs in var- 
ious systems of bosons, including, in addition to atoms, quasiparticles in systems out of 
equilibrium such as excitons and polaritons (for example, [1]). Theory predicts that quantum 
spin excitations in solids, being Bose-quasiparticles, can at certain conditions build the con- 
densate, and magnetic ordering in various sysems can be understood in terms of the BEC of 
these excitations. [2-6] The first experimental observation [7] of magnetic field-induced BEC 
of triplons, that is the spin 5* = 1 quasiparticles, in antiferromagnetic TlCuCls produced a 
diverse research field. [8-17] In this compound, the triplon branches with = —1, 0, 1, are 
separated from the ground state by a relatively small gap A. For this reason, the Zeeman 
interaction in a modest external magnetic field H^xt can close the gap for the = —1 states. 
In contrast to atomic gases, where the total particle number is constant, for triplons it is 
proportional to magnetization M{T,Hcxt) induced by i^cxt- The density of triplons rapidly 
increases with the field, and they undergo the BEC leading to a magnetic ordering. This 
field-induced BEC, which occurs at the scale on the order of few K, has been observed in a 
variety of quantum antiferromagnets.[17] 

The condensate properties crucially depend on interaction of the particles. [18] For the 
atomic BEC at T = the interatomic repulsion can lead to the condensate instability when 
the concentration becomes large enough. [19] Another general feature clearly seen in the 
triplon BEC is the dependence of its physics on the bare dispersion of the quasiparticles 
^k- The non-parabolic bare dispersion of triplons [20] leads to a non-trivial dependence of 
the transition temperature Tc on the concentration p ~ M{T,Hcxt) and, hence, on i^ext- 
The bare dispersion, being itself ifext-independent, determines the interplay of kinetic and 
potential energy of a macroscopic system, and, therefore, plays a crucial role in the BEC 
properties. The effects of the bare dispersion are clearly seen experimentally as the p- 
dependence ~ p'^^'^\ The exponent approaches 2/3 at low concentrations (low 

Tc), [11] as predicted for the parabolic ^k, while at T > 2.5 K, 0(p) is close to 0.5. 

Here we establish theoretical phase diagram of the field-induced triplon BEC based on 
the Hartree-Fock-Bogoliubov (HFB) approximation taking into account also a nonparabolic 
dispersion and determine the fields H^^^ and iJpxt > -^cxt? corresponding to the BEC onset 
and to the instability. A problem in the current theoretical description of the transition at 
i^ext is its predicted discontinuity. We show that this result is an artefact of the conventional 
Hartree-Fock-Popov (HFP) approximation, neglecting the anomalous density terms. When 
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the anomalous density is taken into account, the theory correctly predicts the continuous 
transition. For this reason, the HFB method is more appropriate to study the BEC than 
the HFP one. We find here the stability region of the triplon BEC in the T — Hcy± plane and 
prove that its boundaries strongly depend on the dispersion e^- Results on triplons and on 
cold atoms can be compared to foster the understanding of the similarities and differences 
of their BEC. 

The triplons form a non-ideal Bose gas [4, 7, 20] with contact repulsive interaction de- 
scribed by the Hamiltonian: 

H = j dV{i,\r)lCiiv) + ^-[i^{r)i,{r)f], (1) 

where K, is the kinetic energy operator and g is the coupling constant, and we adopt the 
units /cs = 1, /i = 1, and V = 1 for the crystal volume. Below the critical temperature 
Tc the global gauge symmetry becomes broken as realized by the Bogoliubov shift in the 
field operator: ip{r) = v + Here the condensate order parameter v and ip{r) define the 
density of condensed and uncondensed particles, respectively: po = Pi = (^^(r)^(r)). 
The grand canonical Hamiltonian is: Hg = H — fip, where fi is the chemical potential and 
the total density p = po + pi is uniquely determined by i/cxt- The density p is considered 
as a dimensionless quantity. After the Bogoliubov shift one presents the grand Hamiltonian 
Hg in terms of second quantization operators as Hg = Hq + Hi + H2 + H3 + H^ with: [21] 

Ho = -m + ^, (2) 
H2 = (^k - P + 25'po) 4«k + ^ (cika-k + h.c.^ 



ZJ t t 

^4 = 22^k,p,q^kaJ,aqak+p-q, 



where the prime shows that zero momentum states are excluded. Similarly defined linear 
(Hi) and cubic (H^) terms having zero mean- field approximation (MFA) expectation values 
are omitted. 

Now we implement the HFB approximation [21, 22]: 



4«J,«q«m 4aj,am(apaq) + aqam(akap) + 4«J,(«q«m) - 2pi - cr^, (3) 
where (aj^Op) = ^k.p'^k, (okO-p) = 5k,-pC"k) and are related to the normal pi = y^nk 

k 

and anomalous o"i = ^ (Xk densities. The grand Hamiltonian in this approximation involves 



only zero Hq and second order H2 contributions in Ok, 0^: 

(4) 



Ho = -fipo + 1 pl-2pl- a\ 



Wk^ak + — (akfl-k + h-C; 



where OJ^ = — fi + 2gp and 

X, = 2g{po + ai). (5) 

It follows from (4) that for T > the H2 term is diagonal and hence, the triplon density is 
given by the same formula as in the widely used HFP approximation 

p(T > ^c) = E - 1 ^ ^ e(e.- J)/T _ r 

where /iefr = P- — 2gp. In the BEC regime one performs Bogoliubov transformation 

Ok = Mk^k + ^'k&lk) 4 = ^k^k + ^^k^-k, (7) 

with [6k, &p] = <^k,p5 (^k^^k) ~ {bkb-k) = 0. As a result, the grand Hamiltonian is transformed 
to the Bogoluibov form: 

H = Ho + Yl Ekblb^ + I E(^k - cok), (8) 
k ^ k 

where (&k&k) = nsiE^^, T) = l/[exp(£'k/T) — 1] with the phonon Goldstone mode dispersion 
= u^. — Xf/4. At small momenta, this mode is a collective excitation of the condensate 

carrying spin 5*2 = —1, while at large momenta it becomes the bare triplon mode. 

In accordance with Hugenholtz-Pines theorem [23] at small k the phonon dispersion is 

linear: ^ ck + 0{k'^), where c can be interpreted as the speed of sound. This is achieved 

by setting — Xi/2 = 5k, that is, by: 

p- gpo- 2gpi + gai = 0. (9) 



This choice yields E'k = \l s-k.{sk + ^1) with c = yXi/2m, where m is the triplon effective 
mass. It can be shown [19, 21] that Xi is related to the normal and anomalous self energies 
as En = Xi/2 + /i and Sa = Xi/2, respectively. The quantity Xi plays a special role in our 
analysis: when Xi > 0, the condensate is stable, otherwise it decays due to triplon-triplon 
interaction. [24-26] Below we find Xi in the T — if^xt plane and determine the stable BEC 
region by the condition Xi > 0. 



Using the explicit = J {u^ + Ek)/2Ek and = J (wk — -E'k)/2-E'k, one obtains: 



pi = E(«U) = e(^-^) 



k k k 

where W]^ = nB{E]f_, T) + 1/2. Near the transition, T — )■ Tc the condensate fraction vanishes: 
Po — ^ 0, and Eq.(5) yields Xi = 0. In the triplon BEC, the critical density pc corresponds to 
Peff = 0, i.e. Pc = fi/^g. Therefore, at a given chemical potential /i = gpsHe^t — where 
g is the electron Lande-factor, T^. is determined by: Sk?T-B(^k; ^c) = p/2(?. 

To perform MFA calculations one starts by solving Eqs.(5) and (9) with pi and cti given 
by Eq.(lO). In contrast to the BEC of atomic gases, in the triplon problem, the chemical 
potential is the input parameter, whereas the densities are the output ones. Bearing this in 
mind, we rewrite Eqs.(5) and (9) as: 

Xi = 2p + Ag{ai- pi), p^ = -1 - a^. (11) 

Using dimensional regularization at T = 0, we can find from (10) more explicit expressions 
for the densities 



/d fv 



d^k ^,^1/2" ^^^^ 



where ai{T = 0) = 3pi(T = 0) = V^imXif^^ /Att^, as shown in Ref.[22]. By setting in all 
above formulas cxi = 0, one arrives at the HFP approximation, [7, 20] and particularly 

^[HFP] 

Xr^] = 2p-4,p„ po = ^. (13) 

2^ 

The above Eqs.(ll)-(13) can be applied for any realistic ^k- It is instructive to note that for 
the parabolic dispersion ^k = k^/2m, the BEC can be fully described by only two parameters 
f] = fim^g^ and t = T/Tc with = c{fi/g)'^^'^ /m, c = ir (^-\/2/(73/2(l)) ^ = 2.0867, where 
g3/2iz) is the Bose function. [18] The parameter rj is an analogue of the gas parameter [18] 
of atomic BEC. 

Since the MFA (both HFB and HFP) calculations are based on Eqs.(ll), (12), a question 
about the existence of positive solutions for Xi arises. To analyze qualitatively the existence 
of the physical solutions, we consider T = case. Here, the HFP Eq.(13) is simphfied by 
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FIG. 1: (Color online). Comparison of the HFB (solid) and the HFP (dashed lines) results for 
the triplon density. The HFB approach shows a continuous behavior, which fully agrees with the 
experimental data [11, 12] while the HFP approach leads to the discontinuity. The corresponding 
i^ext are marked near the plots. 

substitution Zhfp = x|^^^y2/i to 1 = ^hfp + '^Z'^Y^y/v/'^'^'^ has physical solutions 
^HFP > for any r/ > 0. This remains valid for all t < 1 at any concentration p. However, 
in the HFB approximation the situation is different: even at t < 1, the physical solutions of 
Eq.(ll) can disappear if rj exceeds a critical value rjc- For example, at T = 0, Eq. (11) for 
Z = Xi/2fi simplifies as 

1^.-^. (14) 

When 7] exceeds rjc = vr^/12, the rhs in Eq.(14) is less than 1 for any Z > 0, therefore, it 
has no positive solutions, and, as a result, Xi acquires an imaginary part. Bearing in mind 
that 7] = fim^g"^ = (/i^^ifcxt — A)m^g'^, one concludes that even at T = 0, if the i^ext is 



strong enough the speed of sound c = \J Zji/m becomes complex and, hence, the BEC is 
unstable. 

To calculate p and Xi one needs the bare e\^. Misguich and Oshikawa [20] demonstrated 
that only with the exact one can explain the overall pc — - dependence. Here we 



apply a similar approach, using a simpler, "relativistic" = + — A, generic 

for systems with gapped spectrum. This choice leads to pc ~ at higher and pc ~ Tj^^'^ at 



lower Ts, respectively. [27, 28] Here the effective exchange J = 2yA/m is chosen to match 
the parabolic and the relativistic ^k at small k. 

In numerical calculations we used parameters by Yamada et al. [11] for TlCuCls: m = 
0.0204 K-\ (i.e. m = 0.261 x 10-^5 g)^ unit cell size 0.79 nm, A = 7.1 K, ^ = 313 K 
and g = 2.06. We neglect a weak renormalization of the model parameters by temperature- 
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FIG. 2: (Color online). Phase diagram for the parabolic (a) and relativistic (b) triplon dispersion 
in the HFB approximation. 

dependent many-body effects, which can slightly shift the stability region boundary, since 
we consider the regime of low T and p. This assumption yields a perfect agreement of theory 
and experiment [20] in a similar range of T and p. We begin with the comparison of the HFB 
and HFP approaches for the density p in a constant i^cxt- Fig.l shows a continuous plot 
of p(T,Hcxt) obtained with the HFB approach, [29] in full agreement with the experiment 
[11, 12] and in contrast to the HFP approximation. In Fig.2 we present the phase diagram 
obtained in the HFB for the parabolic and the relativistic ^k- Solid curves in these figures 
present Tc vs Hcxt obtained from EknB(5k,^c) = pf^d- The dashed lines present the BEC 
stability boundary: there is no Xi > solutions to the gap equations in the regions below 
these lines. Therefore, the HFB approach predicts the existence of a stable (the region 
between solid and dashed lines) and unstable BEC zones (the region below the dashed 
line). As expected, at low T and small i^ext the stability region in Figs. 2(a) and 2(b) is 
the same for both ^k- In general, the relativistic dispersion leads to a narrower stability 
zone than the parabolic one. Note that magnetization measurements on TlCuCls have been 
done for Hcxt between 5.1 and 9 T.[ll, 12] It would be interesting to experimentally study its 
behavior at higher H^xt to explore the instability region. [30] A direct access to the dispersion 
and damping of the phonon-like mode in TlCuCla can be achieved in the inelastic neutron 
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FIG. 3: (Color online). Triplon density as a function of T for relativistic (solid) and parabolic 
(dashed lines) dispersion in the HFB approximation for Hcxt marked near the plots. The plot for 
H = 12.5 T (upper curves) shows two anomalies, one of them caused by the instability. 

scattering measurements. [31] 

Density p as a function of temperature is presented in Fig. 3 for two ifcxt- At relatively 
weak fields, e.g. iJcxt = 7.0 T the magnetization exhibits only one anomaly at T = Tc while 
at stronger one Hcxt = 12.5 T, two anomalies are present. The minimum at the solid line 
at 6.2 K is the onset of the BEC, while the anomaly at T slightly less than 3 K is due to 
the condensate decay. Similar physical behavior can be seen in Fig. 4, which shows the BEC 
fraction po/p x 100%. This fraction is rather large (~ 95% for H^xt = 7.5 T at T = 0) 
and rapidly decreases with increasing the temperature. In both Figs. 3 and 4 the curves 
for Hcxt = 12.5 T start at T ^ 3 K since the BEC is unstable below this T. However, 
Fig. 4 shows that even close to this point the condensate fraction is approximately 70%, and, 
therefore, in the instability zone the condensate can exist for a short time [32] determined 
by the imaginary part of the self energy Xi . This regime will be considered in an extended 
paper. 

In summary, we have theoretically established the phase diagram of the field-induced 
triplon BEC in quantum antiferromagnets in the T — Hcxt plane for a model relevant for 
the TlCuCla compound. Our approach is based on the HFB approximation taking into 
account the anomalous density in the condensate phase. We have shown that (i) at the BEC 
transition the magnetization remains continuous demonstrating a minimum, in agreement 
with the experiment, (ii) in high magnetic fields the condensate becomes unstable due to the 
triplon-triplon repulsion, resulting in interaction of quasiparticles, and found the stability 
boundaries. The non-parabolic dispersion of triplons determined by the crystal structure 
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FIG. 4: (Color online). Condensate density fraction for relativistic (solid line) and parabolic 
(dashed line) in the HFB approximation for //cxt marked near the plots. 

has the crucial effect on the phase diagram by changing the boundaries i^ext ^-iid ifgxt ^"^^ 
making the stability region smaller. 
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